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THE STRUCTURE OF THE M87 JET: A TRANSITION FROM PARABOLIC TO CONICAL STREAMLINES 
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ABSTRACT 

The structure of the M87 jet, from milli-arcsec to arcsecond scales, is extensively investigated, utilizing the 
images taken with the EVN, MERLIN and VLBA. We discover that the jet maintains a parabolic streamline 
over a range in sizescale equal to 10^ times the Schwarzschild radius. The jet then transitions into a conical 
shape further downstream. This suggests that the magnetohydrodynamic jet is initially subjected to the con- 
finement by the external gas which is dominated by the gravitational influence of the supermassive black hole. 
Afterwards the jet is then freely expanding with a conical shape. This geometrical transition indicates that the 
origin of the HST-1 complex may be a consequence of the over-collimation of the jet. Our result suggests 
that when even higher angular resolution is provided by a future submm VLBI experiment, we will be able to 
explore the origin of active galactic nuclei jets. 

Subject headings: galaxies: active — galaxies: jets — galaxies: individual (M 87) 



1. INTRODUCTION 

M 87 is one o f the nearest active galaxies (16.7 Mpc; 
[Jordan et alJl2005h which exhibits relativistic outflows. The 
mass of the central supermassive black hole (SMBH) is mea- 
sured to be 3.2 X 10^ from HST observations of the ion- 
ized gas disk {e.g., Macchett o et aUll997l) . Recent analyses 
of th e stellar kinematics s uggest a larger mass of 6.6 x 10^ 
Mq dGebhardt et al.ll20Tll) . This larger mass gives an appar- 
ent size ~ 8 ^as for the Schwarzschild radius rg. This galaxy 
therefore provides a unique opportunity to study the relativis- 
tic outflow with the highest physical resolution in units of r^. 

Based on the VLA/VLBI observations during the past three 
decades, th e structure of the M 87 jet has been extensively 
examined. lJunor et alTI (119991) found a smooth variation of 
the jet opening angle from 60° at 0.03 pc to smaller than 
10° over ~ 10 pc, indicating that the jet is being collimated 
by the magnetohydrodynamic (MHD) process (Meier etal. 
[200 Ih . One of the most interesting features in the M87 jet 
is the bright kn ot G, lying about 1 " 78 pc) from the core in 
the VLA maps dOwen et al. llT989h . That region has been re- 
sol ved_b^;JheIKTmto a structured complex known as "HST- 
1" (iBirettaet al. 1119991) . It is located around 0.8 - 1.0" (or 
62 — 78 pc) from the core. It consists of bright knots, whose 
(apparent) prop er motions a re superluminal with a range of 
4c — 6c (Birett a et al. Ill999h . A similar velocity pattern has 
been observed at radio frequencies using the VLBA. The com- 
ponent of HST-1 which is furthest upstream (i.e., HST-ld), is 
stationary to within the errors (< 0.25c), and has been identi- 
fied as the origin of the superluminal ejections (Cheung et al. 
120071) . The structure of the jet downstream of HST-1 (1 - 18" 
or 0.1 — 1.5 kpc) can be charac terized by a conic al shape 
with an opening angle of - 6° dOwen et al. II1989I) . HST- 
1 is also a remarkable site for large flaring activity across the 
electromagnetic spectrum from radio thr ough optical to X-ray 
bands (Madrid 2009: Ha rris et al. l '2009). Chandra observa- 
tions show that the flaring event started in 2000, and its bright- 
ness eventually increased more than 50 times. Peak brightness 
occurred in 2005 and an adi abatic compression has been sug- 
gested as a cause of the flare dHarris et al. 120061) . It is thought 



to be one possible sites for TeV 7-ray emission (Harris et al. I 
120091) . Note that all distances are in projection throughout this 
section. 

2. OBSERVATIONS AND DATA REDUCTIONS 
2.1. EVN data 

We conducted the EVN observations towards M 87 on 
7 March 2009 at a wavelength of 18 cm with Cambridge 
(UK), Effelsberg (Germany), Jodrell Bank (UK), Knockin 
(UK), Medicina (Italy), Noto (Italy), Onsala (Sweden), Torun 
(Poland), and Westerbork (Netherlands) stations. Both left 
and right circular polarization data were recorded at each tele- 
scope using 8 channels of 8 MHz bandwidth and 2 bit sam- 
pling. The data were correlated at the JIVE correlator 

A priori amplitude calibration for each station was derived 
from a measurement of the system temperatures during each 
run and the antenna gain. Fringe fitting was performed us- 
ing AIPS. After delay and rate solutions were determined, 
the data were averaged over 12 seconds in each IF and self- 
calibrated using Difmap. 

2.2. MERLIN data 

We conducted the MERLIN observations on 9 March 2007 
at a wavelength of 1 8 cm with Defford, Cambridge, Knockin, 
Darnhall, Jodrell Bank and Tabley stations. Both left and 
right circular polarization data were transferred from each 
telescope with 15 MHz bandwidth. 

Initially, the data were pipelined, and a channel-by-channel 
flux scale derived from 3C 286 observation was applied. After 
this a priori amplitude calibration, we self-calibrated the data 
using Difmap. 

2.3. VLBA data 

We also analyzed archival VLBA data of the M 87 jet 
(BK073). Those observations were carried out on 22 January 
2000 at 15 GHz with aU ten stations of the VLBA and one sta- 
tion of the VLA. 4 IFs, each with an 8-MHz bandwidth, were 
recorded at each telescope. Data reduction was conducted in 
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Fig. 1.- 



(a) MERLIN image of M 87. Contours are plotted at -1, 



1, 2, 

1024 X 3.64 mJy beam^^, which is three-times the residual r.m.s. noise. The 
image is restored with a circular beam of 25 arcsec to emphasize the extended 
structure, (b) EVN image of M 87. Contours are plotted at -1, 1, 2, 1024 
X 1.37 mJy beam~^, which is three-times the residual rm.s. noise. The 
synthesized beam is 17.8 mas X 13.1 mas with the major axis at a position 
angle of 76.1°. (c) VLBA image of M 87. Contours are plotted at -1, 1, 2, .... 
1024 X 0.45 mJy beam~^, which is three-times the residual r.m.s. noise. The 
synthesized beam is 1.22 mas X 0.59 mas with the major axis at a position 
angle of -3.99°. 

the same manner as for the EVN observations. We repro- 
d uced an M 87 image, w hich is consistent with that pubUshed 
in lKovalev eTaP (l2007h . 

3. RESULTS 

We show the images from EVN, MERLIN and VLBA ob- 
servations in Fig. [T] Image qualities are summarized in table 
1. The bright core at the eastern edge of the jet and contin- 
uous jet emission are detected in all images. We define the 
core with a Gaussian model which was fitted to the innermost 
bright region in all the images. The isolated component at the 
distance of 900 mas from the core in the EVN image is the 
HST-1 knot. It corresponds to the first bright component to 
the west of the core in the MERLIN image. This feature was 
not detected in the VLBA image at 15 CjHz, since that data 



were taken before the HST-1 component underwent its flaring 
event. One of the most important achievements in our obser- 
vations is a clear detection of continuous jet emission up to 
500 mas from the core in the EVN image. The continuous jet 
structure between 160 and 500 mas has been suggested in pre- 
vious measurements (Reid et al. 1989; Cheung et al. 2007), 
but here they are robustly detected and resolved with high sig- 
nificance in our EVN observations. 

In order to illustrate the structure of the M 87 jet, we ana- 
lyze the jet emission along a position angle of 338° and fit the 
cross sections of the jet with one Gaussian or multiple Gaus- 
sians for all the images. We define the width of the jet as the 
deconvolved FWHM for the one-Gaussian fitting case. For 
the multiple-Gaussian fitting case, we define the distance be- 
tween the outer edges of the FWHMs as the width of the jet. 
We note that the jet cross sections can be represented mainly 
by two Gaussians when multiple Gaussians were fitted. This 
is consistent with the edge brightened m o del suggested by 
previous observat ions jOwen et al. 1 119891 : iReid et all Il989t 
lJunor et al. |[T999l ). 

We show the radius of the jet as a function of the depro- 
jected distance from the core in units of in Fig. |2l Note 
that the radius as defined perpendicular to the jet does not suf- 
fer from orientation effects. We assume that th e viewing angle 
of the jet is 14° based on a beaming analysis ( I Wang & Zhou I 
120091) . Hereafter, we use the deprojected distance along the 
jet throughout the paper We note that 1 mas corresponds to 
530 Ts in deprojected distance. We measure the radius of the 
jet at the innermost region by using the previous VLBA im- 
age at 43 GHz by .Walker et al. (2008) as well. The jet radius 
increases downstream, as clearly shown in Fig. |2] We found 
that the streamlines of the M87 jet can be divided into two 
different regimes: i) a parabolic shape with a = 1.73 ± 0.05 
(z oc r"-, where r is the radius of the jet emission and z is 
the axial distance from the core), on scales up to ^ 2.5 x 
10^ Tg, and ii) a conical shape with a = 0.96 ±0.1 starting 
at around 2.5 x 10^ r^. It is remarkable that the jet propa- 
gates with a single power-law streamline over more than three 
orders of magnitude in distance. Furthermore, it is notable 
that the two observed power-law streamlines intersect at 2 x 
10^ r^, whic h is close to the Bondi radius of 3.8 x 10^ Vg 
(- 250 pcB (lAUen et al. II2OO6 ). Indeed, this is the first ob- 
servational evidence detecting a transition from parabolic to 
conical streamlines in extragalactic jet systems. 

4. DISCUSSION AND SUMMARY 

4.1. Unconfined Structure: Downstream of HST-1 

We consider a conical streamline (z oc r°, a = 1) for su- 
personic jets. In an adiabatic jet, the internal pressure pjd 
decreases with axial distance z as z^"^^ (F: the ratio of spe- 
cific heats). So, we speculate that the constant expansion of 
the jet radius and a conical structure downstream of HST- 
1 to knot A in the M87 jet, requires the same axial gradi- 
ent fo r the external ISM pressure, pism, as pj^t dOwen et al. I 
119891) . If pism decreases slower than pj^t, then pism > Pjct at 
some distance s o that a recoUimation shock will be triggered 
(ISandersll983h . The self-similar solution of a conical stream- 
line for the magnetized case (with a p urely toroidal fiel d com- 
ponent) requires pism oc z^^^b — 4 (IZakamska et al. 1 12008). 
For a general (non-self-similar) case, 6 > 2 is allowed in 

' The original value is replaced by taking into account of the mass differ- 
ence of SMBH. 
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Fig. 2. — Distribution of the radius of the jet as a function of the deprojected distance from the core in units of r^. We used images obtained by previous 
VLBA measurements at 43 GHz (red circles) and at 15 GHz (orange circles), EVN measurements at 1.6 GHz (green circles) and MERLIN measurements at 1.6 
GHz (blue circles). The jet is described by two different shapes. The solid line indicates a parabolic structure with a power law index a of 1.7, while the dashed 
line indicates a conical structure with a of 1.0. HST-1 is located around 5 X 10^ r^. The black area shows the size of minor axis of the event horizon of the 
spinning black hole with maximum spin. The grey area indicates the size of the major axis of the event horizon of the spinning black hole with maximum spin, 
and corresponds to the size of the event horizon of the Schwarzschild black hole. The dotted line indicates the size of the inner stable circular orbit (ISCO) of the 
accretion disk for the Schwarzschild black hole. 



analytical and numerical models ('Tchekhovs kov et al. II2008I : 
iLyubarsky 2009; Komissarov et al. 2009). X -ray observa- 
tions reveal the ISM properties such as the Bondi radius re ^ 
250 pc and the King core radius 1.4 kpc (see, e.^ . , 

lYoung et al.TIOOllDiMatteo et al. l200llAUen et al. I2006h . 
Thus, the region of conical streamlines in the M87 jet lies be- 
tween re and the marginal radius for the power-law decay 
beyond in the King profile, indicating that the ISM distri- 
bution is essentially uniform. We thus rule out that the struc- 
ture downstream of HST-1 is hydrostatically confined by pism 
in order to conform to a conical streamline. 

In order to possess a conical streamline without any over- 
collimation between knots HST-1 and A, the condition pjct ^ 
Pism should be maintained. Knots HST-1 to A do appear 
to be over-pressur ed with respect to the external pressure 
jOwen et al. II1989I) . However, pjct of the inter-knot regions, 
as es timated by the minim um energy argument for the VLA 
data dSparks et al. Ill996h . appears under-pressured with re- 
spect to Pism, as est i mated by the recent X -ray observations 
(lYoung et al. I I200I iRaffertv et al. I llOO g). One possibility 
is an underestimation of the magnetic field strength when 
the toroidal (azimuthal) components are not considered, cf. 



lOwen et al. I (Il989t) . It has been further suggested that the 
magnetic field energy is at least in equipartition (or even larger 
with a factor of 1 ^ 2) with the ene rgy of the radiating ultra- 
relativistic electrons dStawarz et al.ll2005 ). 

We note that over-pressured knots do appear to have trails of 
statio nary recolhmatio n shocks in purely hydrodynamic jets. 
Falle & Wilson"! ( Il985h performed hydrodynamic simulations 
to apply stationary recoUimation shocks to the observed knots 
at VLA scales under the assumption of the shallow ISM gradi- 
ent (pism oc Stationary features, ho wever, are in conflict 
with the observed large proper motions dBiretta et al. 1119951 
119991) . while the ISM also does not appear to have such a 
gradient. Therefore, we suggest that the highly magnetized 
nature of the jet may be responsible for the conical part of the 
M87jet. 

4.2. Confined Structure: Upstream of HST-1 

We next consider a parabolic streamline (1 < a < 2) 
for supersonic jets. It is shown that the magnetized jet 
can be parabolic in analytical and numerical models where 
the ISM pressure is decreasing a s jPism oc z , 5 = 2 
(iTchekhovskov et al. ,.2008;. Komissarov et al. II2009I) . A self- 
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similar solution, also exists for non-magnetized cases with a 
pure parabolic stre amline (a = 2) under t he same pis,„ depen- 
dence with 6 = 2 (IZakamska et al. ir2008 ). Their solution in- 
dicates pjot < Pism with pressure increasing near the jet edge. 
So far, no analytical or numerical solutions, except a = 2, 
have been derived in the non-magnetized cases. We note that 
even in a hydrodynamic recollimation shock model, the jet 
boundary ("shocked zone" at the interface of the jet and the 
ISM) is expected to be parabolic. Ho wever, the ISM is re- 
quired to be uniform (Pism ^ const.) (iKomissarov & Fallel 
|J997; Nalewajko & Sikora]|2009|). 

In order to have parabolic streamlines, 6 '--^ 2 is required 
in magnetized jets. It is important to examine the distribu- 
tions of Pism- A recent a nalytical model for a giant ADAF 
dNaravan & Fabian1l201 ll) shows h ^ 2.3 from beyond tb 
down to the SMBH, and it is universally derived for slow 
rotation of the ISM (< 30 % of the Keplerian speed) with 
the wider range of the g-y iscosity parameter (0.001 - 0.3) 
jShakura & Sunvaevi Il973h . A spherically sym metric ac- 
cretion with MHD turbulence dShcherbakov 112008 ) produces 
6 ^ 2.1 without a strong dependence on the parameters and 
the equation of state. However, the spatial resolution in cur- 
rent X-ray observations is not high enough to resolve this re- 
gion. 

4.3. HST-1 

The radius of HST-1 appears to be smaller than the size 
estimated by the two power-law streamlines. The maximam 
measured radius of the HST-1 region is 2200 ± 300 and 
mean value is 1300 ± 150 Ts, while the expected radii are 
3000 rs ± 400 for the parabolic case and 4100 ± 200 
for the conical case, respectively. One possibility is that the 
jet is subjected to over-collimation towards HST-1. 

As is seen in Fig. |2] HST-1 is located downstream of the 
Bondi radius (re), where the jet changes from parabolic to 
conical shape. This transition is presumably caused by the dif- 
ferent profiles of the external pressure (changing from steep 
to shallow gradients): The distribution of the external gas fol- 
lows the influence of gravity from the central SMBH inside of 
r-Q, while it can be flattened outside of r-g,- The field strength 
|B| at HST-1 is 0.5 - 2 mG as estimated by the X-ray vari- 
abihty (Perlman et al. Il2003: Harris et al. 2009). Supposing 
|B| ~ 1 mG for our reference, the magnetic pressure is ex- 
pected to be about 100 times larger (i.e., over-pressur ed) than 
die IS M value pism ~ 4.3 x 10"^" dyn cm"^ ( Allen et al. I 
HOO^, rie = 0.17 cm-'"' and kT = 0.8 KeV). 

Another interesting objective is understan ding the produc- 
tion of the superlum inal components (Bir etta et al. 1 119991 : 
ICheung et al. I l2007h . The observed motions are modeled 
by the quad relativist ic shocks in the MHD helical jet 
dNakamura et al. |[201 tf). Recent optical polarimetry requires 
a coherent helical field to explain the distribu tion of polariza- 
tion vectors (perpendicular to the jet) (Perl man et al. Il2003h 
and variations of the radio electric vector positional angles a s 
well as the fractional polarization degree (iChen et al. lEoTlh . 
indicating either a helical distortion to the jet or a s hock prop- 
agating through a helical jet (Perlm an et al. ll20Tl1) . 

An interpretation for the origin of HST- 1 as a hydrodynamic 
recollimation shock has been proposed by Stawarz et al"] 
(12006, hereafter S06) and Bromberg & Levinson (2009, 
hereafter BL09). S06 does not explain the observed gradual 
coUimation, and is inconsistent with our observational result. 
They assume a freely expanding (i.e., a conical) jet which is 
already in the particle-dominated regime upstream of HST-1. 



S06 suggests a shallow external pressure profile pism oc ° 
inside some critical radius of 3" (« 234 pc), as determined by 
a stellar cusp (a br eak in the stellar s urface brightness of the 
M87 host galaxy) (iLaueretal. II1992I) . BL09 also models the 
recollimation shock at HST-1 with the same pressure profile. 

In summary, we propose that a transition from parabolic 
to conical streamlines, occurs via a strong compression, due 
to the change in the slope of the axial profile of pism- This 
condition can be achieved by an imbalance between pjct and 
Pism in the radial direction due to the over-collimation. It 
is remarkable that recent General Relativistic MHD simu- 
lations (iMcKinneyl 120061) reproduce the observed transition 
from parabolic (a ^ 1.7) to conical (a ^ 1.0) streamlines. 
However, the scale of the transition does not match with the 
observations. We suggest that the properties of MHD jets 
should be considered in a realistic ISM environment. 

4.4. Towards the Jet Origin: BH or Disk? 

Recent VLBA observations indicate the VLBI core at 43 
GHz (core43) is located at 20 from the position of the 
SMBH for the viewing angle of 14° (iHada et al. 11201 ll) . With 
the VLBA image, the size of the core43 can be measured to 
be 17 ± 4 at the position angle of 68°, which is approx- 
imately perpendicular to the jet axis. If we assume the core 
is the innermost part of the jet emission, the radius of the jet 
is 9 ± 2 r^. This is in good agreement with the extrapolated 
radius of the jet of 9 ± 2 assuming parabolic streamlines 
in the inner region. If this is the case, the jet propagates with 
a single power-law streamline over more than four orders of 
magnitude of distance. 

The jet opening angle is one of the key parameters in 
the M HD jet acceleration mechanism. iBlandford & Pavng 
(119821 hereafter BP82) suggest that the jet is initiated by flie 
magneto-centrifugal force if the jet inclination angle between 
the penetrating magnetic line of force and accretion disk is 
smaller than 60°, which corresponds to a jet opening angle 
larger than 60°. If we assume that the core43 is the innermost 
part of jet emission, the jet opening angle would be estimated 
to be 46° at 20 r^. Therefore, the core43 would be located 
downstream of the Alfven surface where the magnetic pres- 
sure gradient force plays a role. Thus, future sub-mm VLBI 
projects can access initial acceleration mechanisms by direct 
imaging and measure the deviation from parabolic streamlines 
upstream of the core43. 

Furthermore, as we show in Fig. |2] the parabolic jet can 
be extrapolated to the vicinity of the SMBH. At the "ax- 
ial" distance of z = 3.4 rs from the SMBH, it intersects the 
theoretically-expected, inner stable circular orbit (ISCO) of 
the accretion disk for a Schwarzschild black hole. If the jet 
radius departs from the extrapolation line well beyond 3.4 
Ts, it will indicate that the jet originates in the accretion disk 
(BP82). If the streamline is as extrapolated down to 3.4 r^, 
the jet must originate from the accretion disk (BP82) around 
a spinning black hole and/or the ergos phere of the spinning 
black hole (Bl andford & Znaiek|[T977l hereafter BZ77). If 
we see no departure down to 0.5 from the SMBH, it in- 
dicates that the jet originates from the ergosphere of the spin- 
ning black hole (BZ77). Note that the above discussion for 
prograde spinning black holes, while the upper limit of the 
ISCO is 4.5 Ts for retrograde spinning black holes. Therefore 
tracing back the ridge line and identifying the ejection point 
of the jet is one of the crucial tests to an evaluation of the 
formation mechanism of the jet (BP82 and/or BZ77). More 
fundamentally, it provides an important method for detection 
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TABLE 1 
Qualities of the images 



Observation 


[mas] 


Synthesised beam 
[mas] 


[degree] 


Peak intensity 
[Jy beam~^] 


RMS noise level 
[mJy beam~^] 


Dynamic Range 


EVN 


17.8 


13.1 


76.1° 


1.77 


0.47 


3770 


MERLIN 


25 


25 


0° 


2.00 


1.21 


1650 


VLBA 


1.22 


0.59 


-3.99° 


0.838 


0.15 


5590 



the spin of the SMBlfl and evaluation of the general theory of 
relativity. Such conjectures will be explored by direct imaging 
with future submm VLBI experiments. 
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